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Abstract 

We consider the (l+l)-dimensional J\f = (2, 2) super Yang-Mills theory which 
is obtained by dimensionally reducing M = 1 super Yang-Mills theory in four 
dimension to two dimensions. We do our calculations in the large- A'^c approxi- 
mation using Supersymmetric Discrete Light Cone Quantization. The objective 
is to calculate quantities that might be investigated by researchers using other 
numerical methods. We present a precision study of the low-mass spectrum and 
the stress-energy correlator (T^+(r)T^+(0)). We find that the mass gap of this 
theory closes as the numerical resolution goes to infinity and that the correlator 
in the intermediate r region behaves like r~^'^^. 



1 Introduction 



There is a pressing need to solve quantum field theories in the nonperturbative regime. 
Over the last thirty years a significant amount of progress has been made in this area 
using lattice gauge theory. Many of the most interesting quantities in QCD and elec- 
troweak physics are being calculated to ever increasing accuracy. There remains, how- 
ever, a number of nonperturbative quantities in supersymmetric quantum field theories 
that are interesting in a variety of formal physics contexts but have not been calculated. 
Some of these calculations are now just beginning to be considered. With increased 
interest in the physics of extra dimensions, it is more important than ever to solve 
supersymmetric theories in the nonperturbative regime. 

The progress in putting supersymmetry on a lattice has been rather slow due to 
some critical problems: the lack of translational invariance on a lattice, the notorious 
doubling of fermion states 0], and the breakdown of the Leibniz rule 0. Recently, 
however, some interesting new approaches have shed some light on this issue |P, ^ ^ . 
These approaches make possible the restoration of supersymmetry in a continuum limit 
without fine-tuning of parameters and even without introducing some "sophisticated" 
fermions such as domain- wall ||^ or overlap fermions p, However, these techniques 
seem to be applicable to only some subset of all supersymmetric theories. 

Given this increasing interest and promising new ideas for the realization of super- 
symmetry on a lattice, it is worthwhile to provide some specific, detailed numerical 
results using Supersymmetric Discretized Light Cone Quantization (SDLCQ) |jTD|, ^ 
for the simplest theory for which the new lattice techniques are applicable. SDLCQ is 
a well established tool for calculations of physical quantities in supersymmetric gauge 
theory and has been exploited for many supersymmetric Yang-Mills (SYM) theories. 
The M = (2, 2) theory in 1+1 dimensions in the large- iVc limit is discussed in Ref. |ll2| ; 
however, the published results are primitive compared to what can be obtained today 
because of our greatly improved hardware and software. In this paper we are now able 



to reach a resolution of K = 12, while in Ref. |]T2[ we could reach only K = 5. Here 
we will present new and more detailed results on this theory against which the lattice 
community can compare the results of their new techniques. 

Briefly, the SDLCQ method rests on the ability to produce an exact representation 
of the superalgebra but is otherwise very similar to Discrete Light Cone Quantization 
(DLCQ) [|13|, |1^. In DLCQ we compactify the x~ direction by putting the system on a 
circle with a period of 2L, which discretizes the longitudinal momentum as p'^ = nir/L, 
where n is an integer. The total longitudinal momentum P"*" becomes Kir/L, where 



K is an integer called the harmonic resolution [13|. The positivity of the light-cone 



longitudinal momenta then limits the number of possible Fock states for a given K, 
and, thus, the dimension of Fock space becomes finite, enabling us to do some numerical 
computations. It is assumed that as K approaches infinity, the solutions to this large 
finite problem approach the solutions of the field theory. The difference between DLCQ 
and SDLCQ lies in the choice of discretizing either P~ or Q~ to construct the matrix 
approximation to the eigenvalue problem M^jvl/) = 2P+P"|\I/) = 2P+((5")^/a/2|\1/), 
with = Ktt/L. For more details and additional discussion of SDLCQ, we refer the 
reader to Ref. |Tl|. 

An interesting new result of the calculation we present here is that finite-dimensional 
representations of the SDLCQ with odd and even values of K result in very distinct 
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solutions of the M = (2, 2) SYM theory, which only become identical as K approaches 
infinity. One might initially think that this is a shortcoming of the SDLCQ approach, 
but it turns out to be an advantage because it provides an internal measure of conver- 
gence. 

We will give some numerical results of the low-energy spectrum. There we will see 
that as we go to higher and higher resolutions, we find bound states with lower and 
lower mass. We have seen this behavior in the A/" = (1,1) theory where the lowest mass 
state converges linearly to zero as a function of 1/K. This closing of the mass gap as 
K oo was predicted by Witten [l^ for the Af = (1, 1) and Af = (2, 2) theories. We 
find that in the latter case the convergence is not linear in and, while our results 
are consistent with the mass gap going to zero, they are not conclusive. 

We have also been able to solve analytically for the wave functions of some of the 
pure bosonic massless states, and we will present the exact form of the wave function for 
some cases. We will show that the states must have certain properties to be massless, 
which then enable us to count the number of the states for a given resolution K. In 
addition, we will present the formulae to count a minimum total number of massless 
states. 

Finally, we will look at the two-point correlation function of the stress-energy tensor 
(T++(r)T++(0)). We see the expected 1 /r^-behavior in the UV and IR regions, and, 
interestingly, we find that the correlator behaves as in the intermediate region. 

We know of no predictions for this behavior; however, for J\f = (8, 8) SYM theory there 
is a prediction that this correlator should behave like 1/r^ in the intermediate region. 

The structure of this paper is the following. In Sec. |^ we focus our attention on 
the low-energy states. After giving a quick review of A/" = (2, 2) SYM theory with 
SDLCQ, we give some numerical results for the low-energy states, discuss analytically 
some properties of pure bosonic massless states, and present the formulae to count a 
minimum total number of massless states. We discuss the numerical results for the 
two-point correlation function of the stress-energy tensor in Sec. 0. A summary and 
some additional discussion are given in Sec. |[ 



2 Review of Ar=(2,2) SYM theory 



2.1 Af={2,2) SYM theory and SDLCQ 

Before giving the numerical results, let us quickly review some analytical work on 
Ar=(2,2) SYM theory for the sake of completeness. For more details see Ref. |]T2|. This 
theory is obtained by dimensionally reducing A/'=l SYM theory from four dimensions 
to two dimensions. In light cone gauge, where y4_ = 0, we find for the action 



qLC 

^1+1 



dx^dx tr 



+lid.A^y + gA^J+ + V2gele2MXi, Or] + j [Xj, Xj]' 



where x^ are the light-cone coordinates in two dimensions, the trace is taken over the 
color indices, Xj with J = 1, 2 are the scalar fields and the remnants of the transverse 
components of the four- dimensional gauge field A^, two-component spinor fields 6r and 
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9l are remnants of the right-moving and left-moving projections of the four-component 
spinor in the four- dimensional theory, and g is the coupling constant. We also define 
the current J+ = i[Xi,d-Xi\ + 29]^9r, and use the Pauli matrices /3i = ai, P2 = o"3, 
and €2 = —i<y2- 

After eliminating all the non-dynamical fields using the equations of motion, we 
find for = / dx-T+'' 

= J dx-tT{d-Xid-Xi + ield-OR), (2) 

and 

P- = 9'J dx~tT (^-lj+^J+ - ^[Xi,Xj]' + '-{e2(3i[Xi,eR]f^e2(3j[Xj,en]^ . 

(3) 

The supercharges are found by dimensionally reducing the supercurrent in the four- 
dimensional theory. They are 

g+ = 2^/^ J dx-tT{d.Xj(3j^,u,), (4) 



Q- = 9j dx-ti (^-2^/^J+^e2a^Ur, + 2-'/H[Xi,Xj]{(3i(3je2)a^u,^ 

where a, r/ = 1, 2 and are the components of Or. 

We expand the dynamical fields Xj and Ua in Fourier modes as 



(5) 



y/2n Jo \/2k+ 



1 dk+ 



; - ^ f -^[B^,,{k^)e-^'''^' + Bljk^)e^'^^n, (7) 

where p,q = 1,2, Nc stand for the color indices, and A, B satisfy the usual commu- 
tation relations 

[Ai,,{k+), = 6ij5,r5,s5{k+ - A;'+), (8) 

i?J,,(A;'+)} = 5^p5,r6,s5{k+ - k'+). (9) 

We work in a compactified x~ direction of length 2L and ignore zero modes. With 
periodic boundary conditions we restrict to a discrete set of momenta 

TT /■ vr °° r 

k+ = -k, k = 1,2,3,..., dk+^-Y,, 6{k+ - k'+) ^ -6kk' (10) 
L J L^^^ n 

Relabeling the operator modes ^J^a{k) = A{k+ = ^) and \J^h{k) = B{k~^ = ^), so 
that 

[aipg{k),al^^{k')] = 6ij6pr6gsSkk', {bapq{k),bl^^{k')} = 6ap6pr6gs6kk'- (H) 
the expansion is 

1 00 1 _ J. 

XUx-) = E + aXp^k-^y^^^-'l (12) 
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1 °° 1 + 
In terms of a and b, the supercharges are given by 



(13) 



(14) 



k=l 



and 

Qa 



TT 



TT 



E ^(fcl+fc2),fc3U^2) 



a?7 



(15) 



ki ,fc2,A;3=l 



1 



1 / 



)4mj(^2)&TO(A;3)] 



1 



/C2 



[4im{kl)blmj{k2)aiij{k3) - a\ij{k3)aiim{ki)br,mj{k2)] 

[(llij{k3)bjjim{ki )aimj{k2) - blim{ki)a\^j{k2)aiij{k3)] 



-J-[blij{k3)br^i^{ki)brjmj{k2) + ( /^l ) ^^mi ( ^2 ) ( ^3 ) ] 

Ir.t 



-^[&,iij(^3)fevm(A:i)6^mj(A;2) + &,lim(^l)^',lmj(^2)fcTO(A;3)] 



k3^ 



+ — Hij(k3)brnm(ki)brpnj{k2) + ^im(^l)^mj (^2)&to(A;3)] 



t 



+'^{^2)lj\^-^^J^^[blij{kz)aum{kl)ajrnj^^^ 

+ ["k- ik3)baim{ki)ai^j{k2) + 4im(fci)«Lj(^2)a/ii(fc3)] 



)] 



^^^ [4ij(fc3)ajim(fel)&amj(fe2) + i(^l)ii(^2)ajii(fc3)] j |- 



using the relation ([/?/, 62)077 = (^a7?(e2)7j- 

They satisfy the super algebra conditions for anticommutators involving Q^, 

{g+, g+} = 5a/32\/2P+, {g+, g^} = o. 



(16) 



but do not satisfy the the condition {Q^ , g^ } = Sai3'2V2P . Instead, in SDLCQ we 
find 



{g^, g« } 7^ if a ^ (gr)' = V2P,- + V2P^ = {q. 



(17) 



C(5 s 7- ^ ^ -r y^i J — V ^-t 1 y- w ^± 2 — ^^2 

Although we have different P~ for different g~, we can define a unitary, self-adjoint 
transformation C, such that 



CaiijC — a2ij, CbiijC — —b2ij- 



(18) 



and find that CP{C — P2 ■ Thus the eigenvalues of P~ are the same. We may choose 
either one of the two g~'s, at least for our purposes, and in what follows we will use 
Qi and will suppress the subscript unless it is needed for clarity. 
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The momentum, P"*", is given by 



= ^^Qt f = f E + (19) 

We work with a fixed value of momentum 

P^ = \k, K=1,2,... (20) 

We call K the resolution because larger values of K allow larger values of L while 
leaving the momentum P+ fixed. 

The next thing to note is that there are three Z2 symmetries of Qi. The first one 
is i?i-symmetry, where Ra acts as follows 

aiij ^ a2ij, ha -ba (21) 

The second is ^'-symmetry 

The third is what we call T-symmetry 

O'lij —cbiij-i ba unchanged. (23) 

It is easy to see that under these symmetries Qi is invariant. 
Using the relations, 

RiQiRi = -Qt, TQ^T = -Qa, 

we find 

R^{Qt±Qt)Ri = T{Qt±Qt), T{Qt±Qt)T = -{Qt±Qt). (24) 
Also note that 

{Qt ± Qt, Qt ± Qt) = {Qt, Qt} + {Qt, Qt} + ±HQt, Qt} = ^V2P+. (25) 

We work in a subspace of definite momentum so (Qt^Qt ) must have non zero eigenvec- 
tors. Since and Q~ are fermionic operators we see that a bosonic energy eigenstate 
1^^^)++ which is even under R and T-symmetry, can be transformed into 

= QiiQt ± Qt)\'^B)++, - {Qt + Qt)iQt - Qt)m++ (26) 

which are all degenerate with \'$b)++- One should notice here that we cannot use Qi 
and Q2 at the same time since they do not commute with each other. Thus, including 
the supersymmetry, we have an 8-fold degeneracy. Utilizing the remaining S'-symmetry, 
which does not give us a mass degeneracy, we can divide the mass spectrum into 16 
independent sectors. This significantly reduces the size of the computational problem. 
It will be convenient to refer to bound states of this theory as having S, T, or R even 
or odd parity and to refer to a state as having even or odd resolutions if K is an even 
or odd integer. 
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1 n 


1 1 
1 1 


1 9 


1.308 


4.009 


0.0067 


2.144 


0.0040 


1.415 


0026 


1.040 


0.0018 


0.8188 


12.62 


12.24 


0.6304 


2.514 


0.0060 


1.5999 


0.0038 


1.138 


0.0026 


0.8790 


22.06 


15.04 


1.0813 


2.645 


0.4366 


1.712 


0.0048 


1.212 


0.0026 


0.9312 




15.28 


1.1099 


2.773 


0.6016 


1.729 


0.3515 


1.256 


0.0039 


0.9397 




22.53 


1.5732 


2.807 


0.6308 


1.811 


0.4372 


1.347 


0.3062 


1.0072 



Table 1: The mass squared of the first few lowest massive states in the S'-even 
sector in units of g'^N^/ vr for a series of resolutions K. 



K=A 


5 


6 


7 


8 


9 


10 


11 


12 


1.2009 


3.1876 


0.00674 


1.8427 


0.00440 


1.2687 


0.00302 


0.95786 


0.00217 


1.2009 


3.1887 


0.6402 


1.9305 


0.00538 


1.3266 


0.00317 


0.99795 


0.00218 


12.296 


3.3239 


0.6747 


2.0413 


0.45529 


1.4087 


0.00431 


1.0302 


0.00219 


12.296 


11.489 


0.9900 


2.1415 


0.48010 


1.5107 


0.36858 


1.1036 


0.00356 


19.502 


11.492 


1.0313 


2.3603 


0.55873 


1.5219 


0.38647 


1.1345 


0.32053 



Table 2: Same as Table 1 but for the S'-odd sector. 



2.2 Mass gap 

Tables |l] and |^ show the first few low-mass states. We find anomalously light states 
in the sectors with opposite K and S parity for K larger than 4. Furthermore, the 
number of extremely light states increases by one as we increase K by two. We believe 
that these anomalously light states should be exactly massless states, but for some 
reason there is an impediment preventing SDLCQ from achieving this result. Some of 
the evidence for this comes from a study of the average number of partons {n) in the 
bound states. For example, in the sector with S and K even, for each even integer r 
less than K, there is exactly one bosonic massless state with {n) = r. For K odd we 
do not see massless states of this type, but we do find (n) = r for the anomalously 
light bound states in this sector. This is also the first sign of the distinction between 
representations of the supersymmetry algebra in different symmetry sectors, namely 
those with anomalously light states (with opposite S and K parity) and those without 
anomalously light states (with matching 5* and K parity). 

In our discussion of the mass gap we will not include the anomalously light states 
as part of the massive spectrum for the reason given above. To study the mass gap we 
will look at the lowest massive state in each sector as a function of 1/i^' as shown in 
Fig. |I[ There we also show polynomial fits in all four sectors separately. The fits are 
constrained to go through the origin. The quadratic fits look very good in Fig. |l|(a). 
but Fig. |l](b) required a cubic. The two fits with opposite S and K parity look very 
similar as do the two fits with same S and K parity. In each case we could have fit 
all the points with one curve if we were to include a small oscillatory function in the 
fit. We should note here that oscillatory behavior has been observed before in different 
theories |]TB|, |T^. The explanation given there is that those states which show the 
oscillatory behavior comprise non-interacting two-body states. This, however, does not 
seem applicable in our case since the states in Fig. |l] are the lowest energy states; thus 
there are no lower energy states available to form two-body states. 

The distinct character of the mass gap serves as another piece of evidence that 
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(b) 



Figure 1: Plots of the mass squared in units of g'^Nc/ir for the lowest massive states, 
excluding the anomalously light states, with a polynomial fit constrained to go through 
the origin. The plot in (a) corresponds to the sector where S and K have the same 
parity, and the plot in (b) to the sector where S and K have opposite parity. 



we have two different classes of representations. The data is consistent with the mass 
gap closing to as iiT — >^ cx), especially for the case where S and K have the same 
parity. The odd and even representations approach each other as K increases and we 
hypothesize that they become identical in the continuum limit oi K oo. When we 
present the correlation function in Sec. |^, we will see further evidence for this claim. 



2.3 Massless states 

2.3.1 Pure bosonic massless states 

Let us investigate the properties of pure bosonic massless states in full detail in the 
Nc ^ oo limit. This is done by generalizing the discussion of the bound states in 
SDLCQ for M=il,l) SYM theory, as given in Refs. P, m, to M=i2,2) SYM theory. 



For simplicity, let us consider the states consisting of a fixed n number of partons 
only. A pure bosonic massless state is given by 

1^,0) = N Y: E^fe+...+.n),i^/[^!...^„](?i • ■■QnMaUqi) • • • aljg.)]|0), 

where is the normalization factor, gj = 1, 2, . . . is the unit of the light-cone momentum 
Pi = qiir/L carried by the i-th parton, Ai = 1,2 indicates the flavor index for each 
parton, the sum J2a is the summation over all possible permutations of the flavor 
indices Aj's, / is the wave function, and the trace is taken over the color indices. Note 
that we don't have the symmetry factor coming from the cyclic property of the trace 
in the above notation; one has to put in the symmetry factor by hand if one would like 
it to be in there as we will do so for an example given later in this subsection. In other 
words, Fock states with non-zero symmetry factor are not normalized. 
Due to the cyclic property of the trace, we have 

f[Ai...An]{Ql, ■ ■ ■ ,Qn) = f[A2...AnAi]{Q2, ■ ■ ■ , Qn, Ql) = ■ ■ ■ = /[yl„Ai...A„_i] (g'n, • • • , Qn-l) ■ 
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Since = {Q~y/\/2, all the massless states should vanish upon the action of 
Q~ . Thus, we must have (5~|^',0) = 0. This identity, however, can be simplified 
somewhat for pure bosonic massless states. That is, the terms to consider in Q~ are 
those which annihilate one boson and create one boson and one fermion, and those 
which annihilate two bosons and create one fermion. Both the former and latter class 
of terms in Q~ separately annihilates |^E',0). In the large- iVc limit the former class 
gives, writing /(gi, . . . , g„) = ^gi . . .g„/(gi, • • • , gn), 

° " ) { ^^^"^ '^^^^ f[Z . . (gl > • • • . Qn-l + t, qn) 
"/[a!...A„_i,7](^1' • • • " qn + t) 



+ty 



M' 



-^^^^:P^/[a1...A„_2,7,A„](91' • • • ' ' (In-l + Qn), (27) 



and the latter yields 



-^n—lt-^n k 

(28) 

where M,"f = [{PjPj - PjPi) 

^2\ai3, t is the momentum of the created fermion, and the 
momentum conserving Kronecker's delta is understood implicitly. These 

are the necessary and sufficient conditions for a pure bosonic state to be massless. One 
should notice that the above equations reduce to the corresponding equations found in 
Ref. |ll], |T8[ with {^2)013 = 1, = 1 for all i's, and M"j = 0, as expected. 

In principle, we could find the properties of all kinds of pure bosonic massless states 
using Eqs. ( |27| ) and (pS]). However, we limit ourselves here to the investigation of only 
two special types. To simplify the notation, we omit the superscript (0) from the wave 
function / hereafter. 

The simplest case is where n = K, that is to say, all the partons have one unit of 
momentum n/L and, thus, f = f. In this case Eq. ( PTj) is trivially satisfied since we 
cannot have states with [K + 1) partons. From Eq. (p8|) we get 

= /[yli...A„_2,l,2] - /[Ai...A„_2,2,l], (29) 

where we have omitted (gi, . . . , g„) = (!,...,!). Eq. (^) means, with the help of the 
cyclic property of /, that the wave function is unchanged after moving any fiavor index 
to anj/ location in the list of indices. For instance, we find, writing f[Ai...A„] = [^i ■ ■ ■ ^n], 

[1212] = [1221] = [2121] = [2211] = [2112] = [1122]. 

It is clear that the state with the above six wave functions being the same and all 
others zero satisfies (|^), or equivalently Eqs. (|27|) and (|28|) , the necessary and sufficient 



conditions to be massless. Therefore, writing tr[a]^^(l) . . . aj^^(l)]|0) = Ai . . . An, we 
find the state 

A^[1212](1212 + 1221 + 2121 + 2211 + 2112 + 1122) = iV[1212] (2(1212) + 4(1122)) 
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is massless, where we used the cychc property of /. In terms of the normahzed Fock 
states tr[al^(l) . . . aljl)]|0)/(v^iV,"/2) = Ai . . . An = Ai . . . An/ {^N^;^'^), where s is 
the symmetry factor, we find, after normahzing properly, that 

^(1212) + ^(1122) 

is massless since s for 1212 and 1122 equals two and one, respectively. Indeed we have 
found the very same massless state in our numerical results. 

As we have seen above, there is a one-to-one correspondence between a massless 
state and a given set of flavor indices, which has a fixed number of I's and 2's. This 
means that every time we change the number of I's (or 2's) in the flavor indices, we find 
a new massless state. Since we can have K + 1 such different sets of flavor indices, we 
have K + 1 massless states of this kind. As verification of our argument, we enumerated 
all the massless states for K up to six and found all of them with the correct coefficients 
in our numerical results. 

The next case to consider is where n = K — 1. In this case only one of the partons 
has two units of momentum, so that / = \/2f. However, since all the /'s have the 
same factor of -\/2, we can absorb \/2 into the normalization factor N and practically 
can set f = f ■ We have t = 1 and Qi = 1 with z = 1, . . . , n in Eq. (^) and find, writing 
(gi, . . . , g„) = (1, . . . , 1, 1, 2) = (1, 2) and so on, 

= [Ai . . . An]{2, 1) - [Ai . . . A„](l, 2), (30) 

= [Ai . . . A„_2, A„](l, 2) - [Ai . . . A„_2, An, A„_i](2, 1), (31) 

= [Ai . . . An-2, A.-l](l, 2) + [Ai . . . A„_2, ^n, ^n](2, 1), (32) 

where v4„_i 7^ An in Eqs. (|3T1) and (|3^) . For Eq. (^) we have t = 2, A; = 1, 2 and qi = I 
with i = 1, . . . , — 2, and we get 

= [Ai . . . An-2, A, A](l, 2) - [Ai . . . An-2, A, A](2, 1), (33) 



= [A...A„_2,l,2](l,2) + [Ai...A„_2,l,2](2,l) 

-[Ai . . . An-2, 2, 1](1, 2) - [Ai . . . A„_2, 2, 1](2, 1). (34) 

Apparently, we have five equations for the massless states to satisfy, but it is easy to 



see that Eq. (^) is incorporated into Eq. (^) and that if Eqs. (|30| ) and (^) are true, 
so is Eq. (0) automatically. Hence, the three equations Eqs. (|30|) , (^I]), and (|3^) are 
in fact the equations for massless states to satisfy for n = K — 1. 

In order to see what the three equations allow us to do, let us first write 

[Ai,...,A„](l,2) = [Ai,...,<]. 

That is, let us put a prime on top of an index whose corresponding parton has two units 
of momentum. Then, Eq. ( ^0|) allows us to move the "prime" to any index. Eq. (^l]), 
along with this fact, then also allows us to move the index with a prime to any location 
in the index list. For example, we have 

[112'] = [11'2] = [1'12] = [12'1] = [1'21] = [121'] = [2'11] = [21'1] = [211']. 
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Furthermore, Eq. (p^ allows us to replace 11' by 2'2 (or 22' using Eq. (pOD) as long as 
a minus sign is inserted. Thus, for the above example we get 



[112'] = [11'2] = [1'12] = -[22'2], 

where we have omitted the wave functions related by cyclic permutations. This means 
that the state 

(112' + ll'2 + l'l2-22'2)/2 

is massless. Note that the symmetry factor in this case is equal to one for all the Fock 
states above. 

Since Eqs. (0), ([3T|), and ( p2D relate all the sets of flavor indices with an even/odd 



number of I's to one another, we have only two independent sets of flavor indices: the 
one with even numbers of I's and the other with an odd number. This means that 
there are two massless states of this type. Again we have conflrmed this statement 
numerically for K up to six. 

To summarize, we have found in the large- iVc limit the necessary and sufficient con- 
ditions, Eqs. (pTl) and (^), that pure bosonic massless states are to satisfy. As an appli- 
cation we considered two special cases and found that there are K + 1 massless states of 

the type tr[ai^(l) . . . a\^{l)] and two of the type tr[ai^(l) . . . a\^_^{l)a\^,_^{2)]. Also, 
we gave a way to enumerate all such massless states for a given K. 

2.3.2 Count of massless states 

It is possible to predict a minimum number of massless states by comparing the number 
of states in the different symmetry sectors. Since (<5~)^ takes a state from one symmetry 
sector to another and then back it must have eigenvalues if the dimensionality of the 
intermediate sector is less than that of the original sector. It is possible to create a 
simple recursive formula for the number of states in each sector |TP|. For the case when 



K is prime and odd, the formula is particularly simple. We present the results here but 
refer to the other publication for justification. We define Af,es+{K,n) as the number 
of states in the bosonic sector with an even number of partons and even S symmetry, 
where n indicates how many types of particles we have in a SYM theory, i.e. n = 4 for 
X = (2, 2) SYM. Then 

{K, n)=A fes+ (K, n) = ^ (35) 



Abes- {K,n) = A fes- {K, n 



2 

Abos+{K,n) = Afos+{K, n) = Af,os-{K,n) = Afos-{K, n) 
where 

1 - ,K 



Af{K,n) + Af{K,-n)-W 

Aj{K,n)-Af{K,-n) 



AfiK, n)pri^e = ^((1 + - (1 + n)) (36) 
W={-nK-l) (37) 
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Q goes from bosonic to fermionic and from even to odd. 



W 



n] 



-Af{K,-n)-- (3. 



2 

W 



A fos- {K,n)- Ahes- {K,n) = A^^s- {K, n) - A {K, n) = -Af {K, + y (39) 



The minimum total number of massless states must therefore be 

-AA^iK, -n) = - nf - (1 - n)) = -|(3^ - 3) 



(40) 



For K = 5, this comes to 96 states which is way more than the 8 purely bosonic states 
with 4 or 5 partons that we have found in this section. 



3 Correlation functions 

One of the physical quantities we can calculate nonperturbatively is the two-point 
function of the stress-energy tensor. Previous calculations of this correlator in this 
and other theories can be found in 12^, ^ Ref. |2y] gives results for the theory 
considered here but only for resolutions K up to 6. We can now reach K = 12. 

We will show that there is a distinct behavior for even and odd K in the correlation 
function, just as in the energy spectrum. Then we will argue, by taking a closer look 
at the data, that we have two different classes of representations at finite K, which 
become identical as K ^ oo. 



3.1 Correlation functions in supergravity 

Let us first recall that there is a duality that relates the results for the two-point 
function in A/'=(8,8) SYM theory to the results in string theory [^. The correlation 
function on the string-theory side, which can be calculated with use of the supergravity 
approximation, was presented in pOf, and we will only quote the result here. The 
computation is essentially a generalization of that given in [p3| , The main conclusion 



on the supergravity side was reported in [^. Up to a numerical coefficient of order 
one, which we have suppressed, it was found that 



{oix)om 



(41) 



This result passes the following important consistency test. The SYM theory in two 
dimensions with 16 supercharges has conformal fixed points in both the UV and the IR 
regions, with central charges of order A^^ and iVc, respectively. Therefore, we expect the 
two-point function of the stress-energy tensor to scale like N'^/x^ and N(./x^ in the deep 



UV and IR regions, respectively. According to the analysis of ||2^, we expect to deviate 
from these conformal behaviors and cross over to a regime where the supergravity 
calculation can be trusted. The crossover occurs at x = l/^fyAfV^^ and x = y/l^/gyM- 



At these points, the Nc scaling of (41) and the conformal result match in the sense of 
the correspondence principle . 
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We should note here that this property for the correlation functions is expected only 
for Af={8,8) SYM theory, not for the theory in consideration in this paper. However, 
it would be natural to expect some similarity between A/'=(8,8) and A/'=(2,2) theories. 
Indeed, we will find numerically that (41) is almost true in A/'=(2,2) SYM theory. 



3.2 Correlation functions in SUSY with 4 supercharges 

We wish to compute a general expression of the form F{x^,x^) = {0{x^ ,x^)O{0,0)) 
where O is r++. In DLCQ, where we fix the total momentum in the x~ direction, it is 
more natural to compute the Fourier transform and express the transform in a spectral 
decomposed form pO|, ^ 



F(P_,x+) = ^(T++(P_,x+)T+^ 
E^(0|T++(P_,0)|z)e-^n^"(.|r++(- 



P-,0)|0). 



(42) 



The position-space form of the correlation function is recovered by Fourier transforming 
with respect to P_ = P"*" = Kti/L. We can continue to Euclidean space by taking 
r = a/2x+x~ to be real. The result for the correlator of the stress-energy tensor was 
presented in [PU|, and we only quote the result here: 



P(a;-,x+) = (r++(x)r++(0)) = \-{0\T++iK)\t) 



2 X~^ 



X' 



Mf 



87r2i^3 



(43) 

where x has light cone coordinates x ,x~^, Mi is a mass eigenvalue and K^lx) is the 
modified Bessel function of order 4. In |[T^ we found that the momentum operator 
T++(x) is given by 



tr 



{d-X'Y + -{iu'^d^u'^ - i(9_M°)M") 



La = 1,2, 



(44) 



where X and u are the physical adjoint scalars and fermions, respectively, following 
the notation of [Q. When written in terms of the discretized operators, a and 6, 
(Eqs. (|ll,|l)), we find 



k=l 



-^k{K - k)a\,^{K - k)aUk) +[--k) 6^,(ir - k)bUk) 



K 



|0). (45) 



The matrix element (L/7r)(0|T++(i^)|z) is independent of L and can be substituted 
directly to give an explicit expression for the two-point function. We see immediately 
that the correlator behaves like at small r, for in that limit, it asymptotes to 



X 

x^ 



F{x ,x^) 



N^{2nb + nf) 
47rV4 



1 

K 



(46) 
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On the other hand, the contribution to the correlator from strictly massless states is 
given by 
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i:|-{0|r-(/o|i).,_^„^. (47) 



That is to say, we would expect the correlator to behave like at both small and 
large r, assuming massless states have non-zero matrix elements. 

3.3 Numerical results 

To compute the correlator using Eq. (HI), we approximate the sum over eigenstates by 



a Lanczos p8[ iteration technique, as described in |21|, Only states with positive 
Rai T and S parity contribute to the correlator. The results are shown in Fig. |^, which 
includes a log-log plot of the scaled correlation function 



2 " .2^4 



f = (r++(x)r++(o)) ( — 1 — ^ (4j 



and a plot of (ilog]^Q(/)/c/log]^Q(r) versus log]^Q(r), with r measured in units of \Jti/ g'^N,.- 
Let us discuss the behavior of the correlator at small, large, and intermediate r, sepa- 
rately in the following. 

First, at small r, the graphs of / for different K approach Qa.sK increases. This 
follows Eq. (|46D which gives the form / = log(l — Second, at large r, obviously, 
the behavior is different for odd in Fig. ^(c) and (d), and even in (e) and (f). 
However, the difference gets smaller as K gets bigger, as seen in Fig. |](a). The reason 
for this is as follows. Looking at the detailed information of the computation of the 
correlator, we found that for even K there is exactly one massless state that contributes 
to the correlator, while there is no massless state nor even an anomalously light state 
that makes any contribution for odd K. Instead, it is the lowest massive state that 
contributes the most for odd K. This observation serves as another piece of evidence 
for the claim that we have two distinct classes of representations for odd and even K. 

In the intermediate-r region, for the A/'=(8,8) theory we expected from Eq. (l4ll) 
that the behavior is an 





21 


21 





this behavior. We indicated 
the derivative of the scaled correlator at a value of —1. Our resolution was not high 
enough to see this in the A/'=(8,8) case. Here we find such a flat region, indicating 
that the correlator in fact behaves like i/r~^-'^^ for A/'=(2,2) SYM theory. Also, note 
that the region of flattening around —0.75 extends farther out as K gets bigger, for 
both odd and even K, implying again that the representations appear to agree as K 
goes to infinity. For any fixed value of r the correlators for odd and even K approach 
each other as K increases and the fiat region extends further. This indicates that it is 
only in the region of r where the correlators for even and odd K agree that we have 
sufficient convergence for the results to be meaningful. 

4 Discussion 

To respond to the increasing interest in calculating supersymmetric theories on a lat- 
tice 0, m, ^ , we have presented detailed numerical results for the low-energy spectrum 
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and the two-point correlation function of the stress-energy tensor, using SDLCQ for 
A/'=(2,2) SYM theory in 1 -(- 1 dimensions in the large- A'^c approximation. Our hope 
is that these results will serve as benchmarks for others to compare and check their 
results. 

In addition, we found an important new aspect of the SDLCQ approximation in this 
calculation. There seem to be two distinct classes of representations for Af=[2,2) SYM 
theory, one where S and K have the same parity and one where S and K have opposite 
parity; these representations become identical as K — oo. We found evidence for this 
feature of A/'=(2,2) SYM theory in both the mass spectrum and the correlator. We 
also found that there are some anomalously light states that appear only in the sectors 
where 5* and K have opposite parity. We argued that the anomalously light states 
should be exactly massless, but have acquired a tiny mass because of some impediment 
to having them exactly massless in the SDLCQ approximation. In the calculation of 
the correlator where only positive S parity contribute we found that there is exactly one 
massless state that contributes to the correlator when K has positive parity and that 
no massless state or anomalously light state contributes when K has negative parity. 
The lightest massive state in the sector where K has negative parity does contribute 
to the correlator, but because the mass gap appears to close at infinite resolution this 
state appears to become massless, as expected JlS 



The two-point correlator of the stress-energy tensor was found to show 1 / r^-behavior 
in the UV (small r) and IR (large r, K even) regions as expected. The large r behavior 
for K odd, on the other hand, has an exponential decay. Surprisingly, the correlator 
behaves like at intermediate values of r. In J\f={8,8) SYM theory in 1 + 1 

dimensions, the correlator is expected to behave like in the intermediate region, 
and it is interesting that A/'=(2,2) behaves similarly but with a different exponent. 
We were able to confirm this power law behavior with a fiat region in the derivative 
of the scaled correlator. Previously, in our calculation of the A/'=(8,8) correlator at 
lower resolutions, we were not able to find this flat region. We are hopeful that in the 
near future we may be able to conclusively confirm the behavior in the A/'=(8,8) 
theory. Interestingly, we also note that earlier results seem to indicate the same type 
of odd/even behavior for the Af={8,8) theory. 

Analytically, we investigated the properties of pure bosonic massless states and 
found the necessary and sufficient conditions to determine their wave function. Then 
we explored some special cases to find that there are K + 1 massless states of type 



trl«l(l)4.(l)---«L(l 



where Ai is a flavor index and the number in the parentheses tells how many units of 
momentum each parton carries, and that there are two massless states of the type 

tr[al^(l)al^(l)...al^^._^(2)]|0). 

We also gave the formulae to count a minimum total number of massless states for a 
SYM theory which is demensionally reduced to one spatial and one time dimensions. 

What prevents us from reaching even higher K is obviously the fact that, as one can 
show [ff9|, the total number of basis states grows like ~ {1 + n)^ , where n is the total 



number of particle types and n = A for A/'=(2,2) SYM theory. Our numerical results 
were obtained using one single PC with memory of 4 GB. The problem that we now 
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face is that we do not have enough memory to store all the states in one PC. However, 
as we make use of a cluster of PCs and find ways to split and share the information 
among them, we are able to reach even higher K. This is the direction of our future 
work, with the ultimate goal being to achieve sufficient numerical precision to detect 
the correspondence between M = (8, 8) SYM theory and supergravity conjectured by 
Maldacena . 
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(e) (f) 

Figure 2: Plots of log of the scaled correlation function / as a function of log^Q(r) for 
(a) X = 3,4, . . . 12, (c) K odd, and (e) K even, and plots of d\og^Q{f ) / d\og'^Q{r) as a 
function of logio(r) for (b) X = 3, 4, . . . 12, (d) K odd, and (f) K even. 
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